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PREFACE 


This  paper  was  prepared  by  the  Institute  for  Defense  Analyses  for  the  Advanced 
Projects  Research  Agency  under  a  task  entitled  "Ultra-Wideband  Radar  Technology 
Evaluation." 

This  paper  was  reviewed  by  IDA  staff  members  Dr,  Arthur  Krinitz  and  Dr.  Roger 
Sullivan. 

This  work  represents  a  contribution  to  the  ongoing  effort  directed  toward  the 
improvement  of  our  understanding  of  the  applications  of  ultra-wideband  radar. 
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IV 


SUMMARY 


This  paper  examines  an  important  basic  issue  which  underlies  the  construction  of 
algorithms  used  in  synthetic  aperture  radar  (SAR)  image  reconstruction;  in  particular,  SAR 
image  reconstruction  in  the  context  of  ground-penetrating  radar  (GPR)  applications.  It  is 
currently  the  case  that  the  development  of  GPR  image  reconstruction  algorithms  is  often 
based  on  assumptions  developed  in  the  context  of  free-space  radar  propagation  to  the 
reconstruction  of  images  of  buried  objects  from  the  SAR  data  [1,2].  Although  in  many 
cases  the  assumption  of  free-space  propagation  leads  to  results  which  are  acceptable;  in 
some  important  cases  it  will  result  in  severe  degradation  to  radar  images  of  underground 
objects.  We  conclude  that  the  assumption  of  plane  wave  radar  propagation  gives  rise  to 
poorer  image  quality  as  the  SAR  receiving  antenna  is  brought  into  closer  proximity  with  the 
air-ground  interface.  This  degradation  occurs  because  as  the  interface  is  approached,  a 
surface  wave  component  of  the  total  electromagnetic  field  known  as  the  lateral  wave 
becomes  a  significant  factor  in  determining  target  signal  history  and  radar  response,  and 
this  conqxment  is  entirely  absent  in  the  plane  wave  approximation.  Improvements  in  image 
quality  do  result  upon  modifying  the  assumption  of  free-space  plane  wave  propagation  by 
incorporating  the  effects  of  refraction  on  plane  waves.  However,  additional  improvements 
in  focusing  quality  result  upon  completely  replacing  the  plane  wave  model  with  the  explicit 
solution  dictated  by  the  actual  boundary  conditions  appropriate  to  GPR  problems.  Then  the 
electromagnetic  field  necessarily  includes  the  lateral  wave,  the  analytical  form  of  which  is 
incompatible  with  any  plane  wave  [3].  We  are  therefore  led  to  propose  that,  in  certain 
specific  circumstances,  the  current  reliance  on  plane  waves  as  putative  radar  waveforms  in 
SAR  image  reconstraction  should  be  abandoned  in  favor  or  more  accurate  solutions  to  the 
Maxwell  equations. 


1.  INTRODUCTION 


In  this  paper  we  shall  examine  an  important  basic  issue  which  underlies  the 
construction  of  algorithms  used  in  synthetic  aperture  radar  (S  AR)  image  reconstruction.  In 
particular  we  will  be  concerned  with  S  AR  image  reconstruction  in  the  context  of  ground- 
penetrating  radar  (GPR)  applications.  We  will  study  the  extent  to  which  one  may  apply 
assumptions  developed  in  the  context  of  free-space  radar  propagation  to  the  reconstruction 
of  images  of  buried  objects  ftom  S  AR  data.  It  is  currently  the  case  that  the  development  of 
GPR  image  reconstruction  algorithms  is  often  based  on  such  assumptions  [1,2]. 


2.  ANALYSIS 


We  shall  consider  the  effect  on  SAR  image  reconstruction  of  the  specific  wave¬ 
forms  assumed  to  model  the  actual  electromagnetic  radiation  in  practical  GPR  problems. 
We  shall  take  the  overall  geometry  to  be  modeled  by  two  semi-infinite  half-spaces  (“air” 
and  “ground,”  respectively)  separated  by  a  planar  interface.  We  shall  employ  at  our 
convenience  both  cylindrical  (p,  (p,  z)  and  rectilinear  (x,  y,  z)  coordinate  systems  such  that 
in  each  case  the  z-coordinate  increases  in  value  in  the  downward  direction.  The  ground  will 
be  modeled  as  a  lossy,  electrically  homogeneous  and  isotropic,  non-magnetic  medium.  The 
region  above  the  interface  (air)  will  be  referred  to  as  region  2,  and  the  region  below  the 
interface  (ground)  will  be  referred  to  as  region  1.  The  corresponding  wavenumbers  will  be 
designated  as  k2  (e  R)  and  k\  (e  C),  respectively.  We  will  consider  a  side-looking  SAR 
receiving  antenna  moving  at  a  fixed  altitude  |  z  1  and  displaced  a  fixed  distance  y  from  a 
buried  dipole  located  at  cylindrical  (and  rectilinear)  coordinates  (0,0,d)  (see  Figure  1).^  It 
will  prove  useful  to  define  the  quantities  rq,  ra,  and  ri  as 


(la) 

(lb) 

=  ^p'^  +  (z  +  df  . 

(Ic) 

We  will  also  need  the  function  F(x)  =  ^(\’^i)-C2(x)-iS2(x),  where 
C2(x)-\-iS2(x)  denotes  the  Fresnel  integral  of  the  second  kind  as  a  function  of  x,  and 
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The  analysis  presented  here  is  for  the  one-way  specialization  of  the  complete  radar  i»x>blem  comprising 
a  buried  source  and  an  elevated  receiver.  This  allows  direct  £q)plication  of  prior  work  on  electromagnetic 
propagation  and  can  be  extended  to  the  cmnplete  two-way  radar  propagation  case  [4]. 


We  will  take  the  source  of  the  electromagnetic  field  to  be  a  buried  horizontal  (x- 
directed)  infinitesimal  point  dipole  with  dipole  moment  normalized  to  unity.  The  complete 
solution  to  the  Maxwell  equations  in  region  2  (i.e.,  z  <  0)  for  an  infinitesimal  horizontal 
dipole  with  unit  moment  buried  at  a  depth  d  in  region  1  has  been  worked  out  by  King  and 
Wu  [5]  and  is  given  by 


0.  z)  = 


2rofej 
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Here  (O  is  the  circular  frequency  and  /io  is  the  magnetic  permeability  of  ffee-space. 
The  above  solution  is  valid  everywhere  in  region  2,  subject  only  to  the  restriction 

Ull>  3U2l.2*^ 


In  the  subsequent  analysis  it  will  be  useful  for  us  to  define  the  component  in  the 
direction  of  motion  of  the  SAR  receiving  antenna  of  the  above  solution  as  the  function 

/l(x)  =  £^^(x)  .  (5) 
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Note  that  the  expression  given  for  £]^^(p,0,z)  in  Eq.  (4)  above  differs  from  the  expression  given  in 
Eq.  (5.10.1b)  on  p.  221  of  [5].  There  is  an  incorrect  factor  of  r2  Ip  multiplying  F(P)  in  Eq.  (5.10.1b) 
of  [5]  when  applied  to  [6],  which  should  be  replaced  by  unity,  and  this  is  done  in  Eq.  (4)  above. 

The  inequality  %2|  ‘s  a  consequence  of  the  assumption  that  the  wavenumbers  ki  and  k2  are 

significantly  different  from  each  other,  or  that  |jti|  »  h\- 
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We  are  interested  in  isolating  and  illustrating  the  effect  on  focusing  quality  which  results 
from  successively  improving  upon  the  ffee-space  plane  wave  model  by  first  taking  into 
account  refraction  effects  in  plane  waves  and  then  making  use  instead  of  the  full  King-Wu 
solution.  In  assessing  the  consequences  of  different  choices  of  radar  waveforms  in  SAR 
image  reconstruction  we  shall  thus  introduce  additional  functions  chosen  to  illustrate  the 
successive  improvement  in  focusing  quality: 

/2(x)  s 

[where  r2(z)  is  given  in  Eq.  (Ic)],  which  is  a  plane  wave  in  an  infinite  lossless  medium, 
and 

/3(x)  s  7(01.  W  ,  (7) 


a  refracted  plane  wave  passing  through  lossy  and  lossless  media  in  which  refraction  effects 
and  transmission  effects  at  the  surface  are  accounted  for  through  the  appropriate  refractive 
path  length  and  through  the  transmission  coefficient 


(8) 


In  the  above  expressions 

R  =  ^(p-af  +  ^  , 

and  the  transmission  angles  0i  and  0i  are  given  by 

01  =sin“^|^sin02)  , 

and 

02  =  sin-l(^)  ^ 


(9a) 

(9b) 

(10a) 


(10b) 


where  a  satisfies  the  transcendental  equation  (derived  in  the  appendix) 

q^{a)a^  -2pq^{a)a^  +  -k-lpkld^a-kld^p^  =0  ,  (11) 

with 

q{a)  -  |u|  [Re^i  •  cos(arg  u)  -  Im^i  •  sin(arg  u)]  .  (12) 

The  quantity  u  which  appears  in  Eq.  (12)  is  the  complex  function  of  a  given  in  Eq.  (A4)  of 
the  appendix.  The  solution  of  Eq.  (1 1),  a,  is  the  horizontal  offset  distance  separating  the 


origin  of  coordinates  ftom  that  point  at  which,  in  the  geometrical  optics  approximation,  the 
refracted  ray  from  source  to  receiver  impinges  on  the  air-ground  interface.  Although  Eq. 
(11)  is  very  interesting  to  study  and  appears  to  have  been  derived  here  for  the  first  time  in 
the  litetature,^  we  will  content  ourselves  with  making  use  of  numerical  solutions  for  the 
purposes  of  the  present  analysis. 

Comparing  the  King-Wu  field  given  in  Eq.  (5)  with  the  different  plane  wave 
approximations  given  in  Eqs.  (6)  and  (7),  one  may  note  an  important  distinction:  the  King- 
^^u  solution  includes  a  surface  wave  component  known  as  the  latcvcil  wove,  which  is 
entirely  absent  from  either  of  the  plane  waves.  The  characteristic  features  of  the  lateral  wave 
are  (1)  the  algebraic,  rather  than  exponential,  decay  as  a  function  of  the  radial  coordinate  p 
as  the  field  propagates  along  the  air-ground  interface  away  from  the  buried  dipole  and  (2)  a 
decrease  in  field  strength  as  I  z  I  increases  in  value.  Inspection  of  Eqs.  (2)  through  (4) 
reveals  that  near  the  interface  the  lateral  wave  becomes  an  important  component  of  the  total 
electromagnetic  field.  It  is  therefore  to  be  expected  that,  as  long  as  the  SAR  receiving 
antenna  is  not  too  high,  the  lateral  wave  component  of  the  King-Wu  solution  will  play  an 
important  role  in  determining  propagation  effects  and  therefore  in  determining  radar 
response  as  well.  These  propagation  effects  will  be  markedly  different  from  those 
associated  with  plane  waves.  Although  it  is  the  case  that  the  electromagnetic  field  due  to  a 
point  dipole  in  an  infinite  medium  reduces  in  the  limit  of  infinite  separation  from  the  source 
to  an  infinite  plane  wave,  such  a  plane  wave  can  never  by  itself  include  the  lateral  wave  [3]. 
One  is  naturally  led  to  inquire  as  to  the  extent  to  which  such  differences  will  manifest 
themselves  in  the  quality  of  images  constructed  on  the  assumption  of  plane  wave 
propagation. 

We  will  therefore  study  the  following  functions  to  assess  the  image  quality 
associated  with  the  different  choices  of  radar  waveform  listed  in  Eqs.  (5)  through  (7): 

A,(x)  s  {arg[/i(x)]  -arg[/,(x)]}p^^_^^p^  (14) 


^  In  most  treatments  of  plane  waves  incident  on  plane  boundaries  (see,  e.g.,  [7])  the  point  at  which  (in 
the  geometrical  optics  approximation)  the  geometrical  optics  ray  passes  from  one  medium  to  the  other 
is  taken  to  coincide  with  the  origin  of  coordinates.  In  our  case  this  point  is  in  general  necessarily  a 
function  of  time  dependent  on  the  moving  location  of  the  SAR  antenna,  and  must  be  determined 
accordingly.  Inspection  of  Figure  1  shows  that  the  offset  distance  a  is  determined  by  the  real, 

“geographical”  transmission  angle  §1  rather  than  the  complex,  “formal”  transmission  angle  0i.  The 
relation  between  the  two  transmssion  angles  is  discussed  in  the  appendix. 


and 


(15) 


where  ij  =  The  two  nonvanishing  functions  /t,(x)  and  the  two  nonvanishing 

functions  A,(x)  are  the  magnitude  and  phase  differences,  respectively,  between  the  plane 
wave  approximations  and  the  King-Wu  solution.  In  evaluating  Aj(x)  in  Eq.  (14),  note  that 
the  quantity  inside  the  curly  brackets  has  been  point-wise  adjusted  by  an  appropriate 
multiple  of  2;rto  obtain  the  “phase-unwrapped”  value  of  the  function.  The  three  correlation 
functions  x/*)  1"  Eq.  (15)  (one  auto-correlation  function  and  two  cross-correlation 
functions)  are  the  zero-doppler  cuts  of  the  various  one-way  radar  ambiguity  functions  [8] 
for  the  different  possible  waveforms. 


3.  NUMERICAL  COMPUTATIONS 


We  shall  numerically  study  the  functions  /i,-,  A/,  and  “jQ  for  four  distinct,  representative 
GPR  scenarios.  If  we  denote  the  linear  frequency  by  f,  these  scenarios  are: 

Case  la:  /=  600  MHz,  y=|zl=  3m 

Case  lb:  /=  600MHz,  y=|zl=  30m 

Case 2a:  /=  5  GHz,  y=|2|=  3m 

Case  2b:  /=  5  GHz,  y=lz|=  30  m. 

The  lower  altitude  scenarios  are  intended  to  model  GPR  applications  which  involve  ground- 
based  receiving  platforms,  while  the  higher  altitude  scenarios  are  intended  to  model  GPR 
applications  which  involve  airborne  receiving  platforms.  For  each  scenario,  we  shall  take  for 
region  1  the  conductivity  <Ti  and  the  real  part  of  the  electrical  permitivitty  e  ir  to  be  given  by  ai 
=  0.01  S/m  and  eir  =  9,  respectively,  and  the  corresponding  quantities  for  region  2  shall  be 
given  by  the  respective  free-space  values.  As  noted  in  Section  2,  we  shall  ignore  magnetization 
effects.  Also,  as  noted  in  Section  2,  we  shall  consider  a  one-way  specialization  of  the  complete 
radar  problem  by  taking  the  total  electromagnetic  field  to  be  given  by  the  electromagnetic  waves 
emanating  from  a  buried  dipole  target;  for  all  four  scenarios  we  shall  take  the  dipole  to  be  buried 
at  a  depth  d  =  0.5  m.  For  both  frequencies  at  each  altitude,  graphs  of  the  various  functions  were 
computed  over  a  synthetic  aperture  of  size  comparable  to  the  altitude.  In  all  four  cases  the  data 
are  plotted  over  a  cross-range  “window”  which  is  smaller  than  the  actual  synthetic  aperture:  of 
10  m  width  and  3  m  width,  for  3  m  altitude  data  at  600  MHz  and  5  GHz,  respectively,  and  of 
20  m  width  and  10  m  width  for  30  m  altitude  data  at  600  MHz  and  5  GHz,  respectively.  The 
plots  are  shown  in  Figures  2  through  13;  their  implications  are  discussed  in  the  following 
section. 


power  (dB) 


Field  Phase  Comparisons  -  f=600  MHz,  y=zp=3m.  d=0.5  m 
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phase  defect  (radians) 


Field  Phase  Comparisons  -  f=5  GHz,  y=2p=3m,  d-0.5  m 


Reid  Phase  Comparisons  -  f=5  GHz,  y=zp=30m,  d=0.5  m 
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Field  Magnitude  Comparisons 

King-Wu  solution  vs.  refracted  and  free-space  plane  waves;  f«600MHz;  y«z*“3m;  d=0.5m 
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Figure  10.  Field  Magnitude  Differences;  Irl  s  3  m;  f  s  600  MHz 


Field  Magnitude  Comparisons 

l-Wu  solution  vs.  refracted  and  free-space  plane  waves;  f«600MHz;  y«z*“30m;  d«0.5m 


cross-range  (meters ) 

Figure  11.  Field  Magnitude  Differences;  Izl  =  30  m;  f  s  600  MHz 


Field  Magnitude  Comparisons 

King-Wu  solution  vs.  refracted  and  free-apace  plane  waves;  £=5GH2;  y«2*«3m;  d“0.5m 
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Figure  12.  Field  Magnitude  Differences;  I  zl  s  3  m;  f  =  5  GHz 
Field  Magnitude  Comparisons 


King-Wu  solution  vs.  refracted  and  free-space  plane  waves;  f^SGHz;  y«z*«30m;  d«0.5m 


Figure  13.  Field  Magnitude  Differences;  Izl  s  30  m;  /  s  5  GHz 


4.  FINDINGS 


The  numerical  results  allow  a  comparative  examination  of  the  functions  A,-,  and 
Xj,  where  the  phase  and  magnitude  differences  are  interpreted  as  predictors,  and  the 
correlation  functions  as  direct  measures,  of  image  reconstruction  quality  [9,  10,  11]. 
Inspection  of  the  correlation  functions  in  Figures  2  and  4  reveals  that  at  frequencies  of 
600  MHz  and  5  GHz,  the  3  m  altitude  plane  wave  images  are  considerably  degraded 
compared  to  the  images  correlated  with  the  King-Wu  solution.  This  is  especially  true  in  the 
case  of  images  correlated  with  the  free-space  plane  wave  (FSPW)  given  in  Eq.  (6), 
whereas  use  of  the  refracted  plane  wave  approximation  (RPW)  of  Eq.  (7)  affords  an 
obvious  improvement  over  FSPW  in  terms  of  both  central  peak  sharpness  and  relative 
height  of  the  nearest  sidelobes.  However,  Figures  3  and  5  reveal  a  significant  difference  in 
the  distribution  of  power  with  respect  to  cross-range  compared  to  the  3  m  altitude  results. 
In  sharp  contrast  to  the  low  altitude  cases,  at  frequencies  of  600  MHz  and  5  GHz  the  30  m 
altitude  FSPW  and  RPW  images  nearly  approximate  the  images  correlated  with  the  King- 
Wu  solution.  Compared  to  the  King-Wu  images,  there  is  a  loss  in  image  intensity  for  both 
FSPW  and  RPW.  However,  at  both  frequencies  the  figure  of  merit  afforded  by  the 
combination  of  central  peak  sharpness  and  relative  height  of  nearest  sidelobes  is  greatly 
improved  for  the  plane  wave  approximations  by  receiving  data  at  a  higher  altitude. 

It  has  been  known  for  many  years  that  the  phase  history  of  a  target  is  generally 
considerably  more  important  than  the  spectral  magnitude  in  the  Fourier  representation  of 
signals  used  to  reconstruct  radar  images  [12].  This  expectation  is  borne  out  by  our  results. 
Figures  6  through  9  show  the  comparative  phase  differences  between  the  plane  wave 
approximations  and  the  King-Wu  solution  for  the  scenarios  considered.  What  is  of  interest 
here  is  the  dependence  on  cross-range  of  the  phase  differences,  since  a  constant  phase 
difference  would  be  unimportant.  The  slope  of  the  phase  defect  furnishes  a  relevant  figure 
of  merit  in  these  graphs.  For  each  of  the  four  scenarios  considered,  the  phase  difference 
function  predicts  the  trend  indicated  by  the  correlation  functions.  For  the  low  altitude  cases 
at  both  frequencies,  the  FSPW  phase  defects  vary  through  as  much  as  2  radians,  while  the 
RPW  phase  defects  are  reduced  at  their  maxima  at  each  frequency  by  a  factor  of  roughly 
one-third,  consistent  with  a  moderate  improvement  in  image  quality.  However,  when  the 
receiver  is  elevated  to  30  m  a  drastic  reduction  in  the  phase  defect  is  observed.  At  600  MHz 


the  maximiim  phase  defect  for  both  FSPW  and  RPW  is  reduced  by  a  factor  of  roughly  60, 
and  at  5  GHz  the  maximum  phase  defect  for  both  plane  waves  is  reduced  by  a  factor  of 
roughly  100,  the  greater  improvement  in  the  latter  case  resulting  fiom  the  greater  number  of 
wavelengths  separating  receiver  and  source.  The  differences  in  field  magnitude  shown  in 
Figures  10  through  13,  although  conforming  to  the  trend  established  by  the  correlation 
functions,  are  as  expected  much  poorer  predictors  of  image  quality  than  the  phase 
difference  functions. 


5.  CONCLUSIONS 


In  sum,  we  have  shown  explicidy  that,  when  the  receiving  antenna  is  very  close  to 
the  ground  above  a  buried  source,  the  use  of  plane  wave  approximations  to  correlate 
signals  will  lead  to  severely  degraded  images  compared  to  images  correlated  with  accurate 
analytical  solutions  to  the  Maxwell  equations  [Eqs.  (2)  through  (4)].  In  contrast,  for  many 
other  cases  of  interest,  when  the  receiving  antenna  is  sufficently  high  (over  30  m  altitude), 
the  quality  of  images  correlated  with  a  plane  wave  approximation  is  quite  good  apart  from 
some  loss  in  image  intensity. 
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APPENDIX 

DERIVATION  OF  EQUATION  FOR  REFRACTED 
PLANE  WAVE  OFFSET  DISTANCE 


If  ©1  and  ^  are  the  formal  transmission  angles,  then  according  to  Snell's  Law 


^isin^i  =^2sin^2  •  (Al) 

•  In  the  problems  considered  in  this  paper  one  has  k\yB\e  C.  However,  to  obtain  the  path 
length  and  transmission  coefficient  for  refracted  plane  waves,  we  are  interested  instead  in 

e  R,  the  real,  “geographical”  angle  which  specifies  through  the  offset  distance  a  the 
location  at  which,  in  the  geometrical  optics  limit,  the  refracted  plane  wave  impinges  on  the 

•  air-ground  interface.  d\  is  determined  by  (cf.  [7],  p.  502) 


sin^j  = 


k2S,m62 


^^+k2Sin^62 


(A2) 


which  is  the  modified”  Snell's  Law  appropriate  to  the  case  of  plane  wave  incidence  at  a 
boundary  separating  lossy  and  lossless  media.  Here  q  is  given  by 


q(a)  =  InfRe^j  •  cos(arg  u)  -  Imki '  sin(arg  n)]  , 

(A3) 

• 

where  i/  is  a  complex  function  of  a  (through  its  dependence  on  &2)  given  by 

u  =  ^-f^a^-/3^-2afiijsin^e2  . 

(A4) 

and 

• 

a  =  ^  , 

(A5) 

l^lP 

(A6) 

• 

Also,  from  Figure  1  we  see  that 

sin02=f  = - 2-^75-  . 

(A7) 

A-2 


We  require  that  the  constraints  in  Eqs.  (A2)  and  (A7)  both  be  satisfied.  We  solve  the 
constraints  by  substituting  for  62  from  Eq.  (10b),  which  leads  to  a  condition  where  a 
satisEes  the  following  transcendental  equation 

q^{a)a^  -  2pq^{a)a^  +  +  Ipk^d^a  -  k^d^p^  =  0  ,  (A8) 


which  is  the  same  as  Eq.  (1 1)  in  the  text. 
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